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pression of excited state contributions to the Bg t mixing parameter and speculate about possible 
explanations. 
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1. Introduction 

Currently, theoretical determinations of the mixing parameters Bb ci (q = d,s) are becoming 
more and more urgent in relation to the Unitarity Triangle analysis. The B-parameters are defined 
as the relative deviations from the Vacuum Saturation Approximation (VSA) of the matrix elements 
of AB = 2 four-fermion operators between B-meson states, i.e. 

(B®\ & t ^ B=2 \B®) 

B Bq = % " 2 q , ^ 2 =[W/z(l-7 5 )^][»(l-75)^] • (l.D 

3JB q m B q 

They encode the low-energy information related to particle-antiparticle oscillations and are for- 
mally accessible to lattice QCD simulations. Nevertheless, a direct computation of Bb q is ham- 
pered by the presence of the large value of the &-quark mass, which imposes the adoption of tiny 
lattice spacings (a <C 1/ (5GeV)) in order to avoid large lattice artefacts. A possible way out is to 
expand the B-parameters in Heavy Quark Effective Theory (HQET), i.e. in inverse powers of the 
ft-quark mass. The leading contribution, also known as the static approximation, is expected not 
to be far from the relativistic value, as previous lattice results have shown. Even so, the naive lat- 
tice discretization of the effective four-fermion operators of the static theory, based on Wilson-type 
light fermions, is affected by a non-trivial renormalization mixing, due to the explicit breaking of 
chiral symmetry, which pushes the numerics up to the edge of our current technology. Although 
Ginsparg-Wilson fermions appear as the natural discretization to study left-left four-quark opera- 
tors, we follow a computationally cheaper approach, based on twisted mass QCD (tmQCD) [0], 
which allows for purely multiplicative renormalization at the same computational cost as with Wil- 
son quarks. 



2. Computational strategy 



Our starting point is the equation relating the left-left operator of the fully relativistic theory 
to the four-fermion operators of HQET, 



\m b ) = Ci(m b ,n)Q l (n) + C 2 (rn b ,n)Q 2 (n) + 



(-)■ 

\m b J 



Qi 



ostat 

VV+AA 



ystat 
SS+PP 



(YhYq) (YhYq) + (VhYsYq) (VhYsYq) ■ 



(2.1) 
(2.2) 
(2.3) 



Eq. ( J2. 1[ ) has to be understood as a scheme dependent perturbative matching between two renor- 
malizable field theories. The coefficients C,, known at NLO in the MS/NDR scheme [0], provide 
the RG evolution from the defining scale m b of the effective theory down to a scale \i 1 GeV. 
Although a natural hierarchy \i <m b has to be assumed in the matching equation, it should be 
observed that, due to the renormalizability of the static theory, the four-fermion operators Q 12 are 
perfectly defined at any scale. In particular, they can be perturbatively evolved up to the RGI point 
through the appropriate 2x2 static anomalous dimension matrix, i.e. 



QT 
QT 



c(n) 



<2i(A0 
(kin) 



(2.4) 
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where c(jti) will be defined later. 

The advantage of introducing RGI operators is twofold. On the one hand, they are truly non- 
perturbative quantities, free of systematic uncertainties related to perturbative truncations. On the 
other, they are regularization independent. As such they can be linked to any specific lattice regular- 
ization, to be chosen on the basis of computational convenience. A simplification of the renormal- 
ization pattern is achieved if we perform a change of basis, i.e. we introduce the primed operators 



QT 

QT + 4Qf 



1 
1 4 



QT 
QT 



QT 
QT 



(2.5) 



This redefinition becomes particularly advantageous on the lattice if the relativistic degrees of 
freedom are discretized according to tmQCD at full twist, i.e. with twist angle a = n/2 [jj]]. 
In particular, from now on we consider the specific case of the B 5 -meson, for which we assume a 
fermion content made of a static quark plus a twisted strange quark belonging to a fully twisted 
(c,s)-doublet. Lighter degrees of freedom, i.e. the up and down quarks, do not need to be further 
specified, since they do not enter the valence sector 1 . In tmQCD the operators Q\ 2 are mapped 
onto their odd parity counterparts £2[ 2 , which renormalize purely multiplicatively, as proved in 
[S]. In other words, with some abuse of notation 



(Gf GI ) 
(Gf GI ) 



limZ( 

a->0 



1 RGI C?oO)) (£[ (#)) tmQCD 



limZo 

a->0 



2.RGI (go(a)) (^2("))taQCD 



(2.6) 



where 



'VA+AV 



= {WhYtiYq) {WhY^Yq) + (WhYfiYsYq) {YhYtiYq) , (2.7) 
■ = {YhYixYq) [YhYixYsYq) + {YhYixYsYq) (VhYnYq) + 

^{VhYsYq) (YhYq) + (YhYq) (YhYsYq)} ■ (2.8) 



The RGI renormalization constants Z' k RGI (k = 1,2) have been recently obtained in the quenched 
approximation [Q] through finite size techniques based on the Schrodinger functional Since 
the latter allows for the adoption of mass independent schemes, the computation of Z' k RGI has been 
performed with standard (untwisted) Wilson fermions. A preliminary study of the non-perturbative 
renormalization for Nf = 2 has been also presented at this conference [0]. 



3. Non-perturbative renormalization in the Schrodinger functional 

In order to study the renormalization of the four-quark operators, we consider a theory with a 
light quark sector consisting of two massless 0(a) improved Wilson-type quarks (yi, Y2) entering 
the four-quark operators, plus a third light spectator quark 1//3 , regularized in the same way, whose 
role will be clarified in a moment. Suitable renormalization conditions can be specified in terms of 

1 This freedom allows to extend the present strategy to Nf = 2 with any kind of dynamical sea, without incurring in 
mixed action issues, such as the adoption of different lattice regularizations for valence and sea quarks. 
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SF correlators made of bilinear boundary source operators 2. SlS2 , 5^ (lying resp. on the two time 
boundaries xq = and xo = T) 

z slS2 [T] =« 6 £C, I (x)rC, 2 (y) , 4 lS2 [r] = a 6 £C,' 1 (x)rc;(y) , (3.1) 

*J x,y 

and the four-fermion operators J2[ 2 . Due to the flavour and parity structure of J2[ 2 , zero-momentum 
correlators need at least three bilinear boundary sources. Two bilinears are placed at xq = and the 
third one at xo = T. Their product gives rise to a generalized source 

w[r l ,r 2 ,r 3 }=j: Vl [r l }L 23 [r 2 }L , 3 - h [r 3 } , (3.2) 

which is parity-odd under five different choices of the Dirac matrices Ti, F 2 and F 3 , i.e. 

^ = W[y 5 ,y 5 ,y 5 ] , sr™ = ~ £ e klm W[y k ,y,,y m ] , 

k,l,m=\ 

^= l -j^W[y 5 ,y k ,y k ], ^ = I £ W[y k ,y 5 ,y k ] , «5^ 5 ) = \ Y.W[y k ,y k) y 5 ] . (3.3) 
-'^i ■ 5 jt=i ■ 5 *=i 

All of the above sources have the same quantum numbers as J2[ 2 and can be used as probes within 
the correlators 

^\x Q )=L-\^ k (x)^) . (3.4) 

Nevertheless, their renormalization is non-trivial and requires the introduction of multiplicative 
renormalization constants to absorb the additional logarithmic divergences of the boundary fields 
from Eq. (3.4). To avoid this, we introduce some boundary-to-boundary correlators 

2Z 6 \^lh 
1 

■ilfel^iW), (3-7) 
0L k=\ 



fi = -^(KM^M) • (3.5) 

fx = "^(^M^iM) , (3-6) 



1 3 



and use them in the ratios 



in such a way that the additional renormalization factors of the boundary sources in Eq. (3.4) drop 
out. The parameter a in the exponent of f[ ! and k l ( can be freely chosen without changing the 
flavour content of the denominator and, in what follows, it will take values a = 0, 1 /2. 

Renormalization conditions, formulated in terms of the ratios h^) a (xo), read 

3gW = l/L) *g,(r/2) =4i(772)L=o , (3.9) 
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where T = L, no background field is introduced and the SF 0-angle [|7|] is set to 6 = 0.5. In our 
simulations we adopt four different lattice discretizations of the heavy quark action, i.e. the stan- 
dard Eichten-Hill one [^] and its statistically improved versions where the naive parallel transporter 
is replaced by a smeared APE, HYP1 or HYP2 gauge link However, in this talk we only report 
on results with the HYP2 action, i.e. the one with the best signal-to-noise ratio. Out of the plethora 
of renormalization schemes that can be defined by Eq. Q, we choose our preferred ones to be 
(s, a) = (1,0) for k = 1 and (s, a) = (3,0) for k = 2 (c/ref. [Q] for further details), thus eliminating 
the indices s and a from the notation. 



4. Renormalization group running 

The formal solution of the Callan-Symanzik equation relates the scheme-dependent RG run- 
ning operator ^(/i) to the renormalization group invariant one (=S[) RGI 



(^)rgi — ^(aO 



POO 

An 




JSMcM, (4-1) 



where g(n) is the scheme and scale-dependent renormalized coupling. Our goal is to compute 
c' k {pi) non-perturbatively. In practice the strategy we follow is to split perturbative and non- 
perturbative contributions at a high renormalization scale pL vU 

(^)rGI = 4(Mpt)^(Atpt,Mhad)^[(Mhad) , (4-2) 

where U k (/J, pt , /ihad ) = c' k (jUhad ) / c' k (jUp t ) represents the evolution of the renormalized operators (ju ) 
from the low-energy hadronic scale /^ad to the high-energy perturbative scale /x pt S> jUhad- Our first 
task has been to compute it non-perturbatively. Since it is difficult to accommodate scales which 
differ by orders of magnitude in a single lattice calculation, it is useful to factorize the evolution 
and adopt a recursive approach. Accordingly, we introduce the so-called step-scaling functions 
(SSFs) o k and a, which describe the change in the renormalization constants and the gauge cou- 
pling respectively, when the energy scale /J. is decreased by a factor of two, 

a(u) =g 2 {n/2) , u = g 2 (n), 

all/2) m=0 

a k {u) = U k (n,n/2y = lim /y P ' '> ^ = limL t ( Bj a M ) , (4.3) 



and go denotes the bare coupling. Having computed the SSFs for a sequence of couplings i = 
0, 1,2, ... 1, we can construct the non-perturbative evolution f/^(2"jUhad, Mhad) from the product 
of SSFs 

£/K2"jUhad I Mhad) = |n (T *( I< 0| . "/=f 2 (2 (!+1 Vhad) • (4.4) 

In the present computation /ihad is taken to be a few hundreds of MeV and we have chosen n = 8, 
so that we could trace the evolution non-perturbatively over three orders of magnitude. In this way 
/i pt = 2" i Uhad is large enough to allow for a perturbative evaluation of c' k {}l pi ) with the operator 



5 



Preliminary results o/Bb s in the static limit from tmQCD 



Filippo Palombi and Mauro Papinutto 



anomalous dimension approximated at NLO [Q]) and the j3 -function at NNLO [11]. The relation 
between the RGI operators and the bare lattice ones defines the total RGI renormalization factor 



(^)rgI =^,RGlfeo)^(«) • 

A comparison between Eq. ( |4~2| ) and Eq. (4^5) leads to 



/t.RGlUW 



: C'k (j"pt ) U' k ( JUpt , Mhad ) 3% (g Q , ajUhad ) 



(4.5) 



(4.6) 



The factor ^'(go^Mhad) must be determined for each operator in a lattice simulation at fixed /^had 
for a range of bare couplings, using suitable renormalization conditions. In our simulations we have 
Mhad = 1/ (2L max ) 270 MeV where L max is fixed through the condition gl P (l/L max ) = 3.480. This 
corresponds to having L max /ro = 0.718(16) (ro = 0.5 fm) JT0|]. The sequence of couplings m; = 
gg F (2~'L max ) is obtained by solving the recursion relation uq = 3.480, c(m/ + i) = w/, / = 0, 1, . .. . 



5. Continuum extrapolation of the step scaling functions 

The lattice SSFs must be extrapolated to the continuum limit (i.e. vanishing a/L) at fixed 
gauge coupling in order to obtain their continuum counterparts at- Since the four-fermion operators 
have not been improved, we expect the dominant discretization effects to be 0(a); therefore our 
data should exhibit a linear behaviour in a/L. Accordingly, we have fitted to the ansatz 

Zk{u,a/L) = a k (u)+p(u)(a/L). (5.1) 

Fits have been performed using either four values of the lattice spacing, i.e. L/a = 6,8,12,16 
or, alternatively, without taking into account the coarsest data L/a = 6, which may be subject to 
higher-order lattice artefacts. The results from three- and four-point fits are always compatible 
within one standard deviation for all operators and schemes, save for a few exceptions in which 
the agreement drops at the level of 1.5 standard deviations only. We have therefore decided to 
choose the three-point based linear extrapolations to extract our final estimates of o k - The resulting 
continuum limit extrapolations for our preferred renormalization schemes are illustrated in Fig. [[]. 
The maximal statistical uncertainty for a* is ~ 1.5% and is found at the largest value of u when 
discarding data at L/a = 6. The values of p obtained in the fits of £l\ are always compatible with 
zero within the statistical uncertainty, while in the case of J2' 2 they are not compatible with zero for 
u>2, thus signalling a stronger dependence upon the cut-off. 



6. Non-perturbative RG running in the continuum limit 

In order to compute the RG running of the operators in the continuum limit as described 
in section ^ we need to fit the results for Ok(u) to some functional form. We follow the same 
procedure as for the renormalized quark mass Jl2|], i.e. we adopt the polynomial ansatz 

M 

<j k (u) = 1 + ^ s m u m , (6.1) 

m=l 
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Figure 1: Continuum limit extrapolation of the SSFs Ei (with s = 1, a = 0) and £2 (with s = 3, a = 0) at 
various SF renormalized couplings computed using the HYP2 lattice discretisations of the static action. The 
SF coupling u increases from top-left to bottom-right. 



with M = 2, 3,4 and s\ always (S2 possibly) set to its perturbative value 



*i = y; (0) ln2, 



s 2 = 7 ; (1) ln2 + 



1 / /(0k 2 , , /(()) 



(ln2) 2 



(6.2) 



It is worth mentioning that if S2 is fitted as a free parameter, it turns out to lie in the ballpark of 
perturbation theory. The RG running factor c[(/ihad) = c[(2"/i had )6 r /(2"/^h a di AVd)> which is now 
a function of the fit parameters only, can be obtained with a complete control of the systematic 
effects. We have indeed checked that its value is fairly insensitive to the fit ansatz and to whether 
S2 is set to its perturbative value or not. We choose to quote as our final results those obtained with 
M = 3, Ji fixed by perturbation theory and s%, S3 kept as free parameters. 

In practice, due to constrains imposed by Heavy Quark Spin Symmetry, the number of inde- 
pendent SF schemes for £!' k is downgraded to four for k = 1 and to eight for k = 2. These lead to 
total RGI renormalization factors which are scheme independent apart from O(a) lattice artefacts. 

The main criterion to define suitable schemes amounts to checking that the systematic un- 



certainty related to truncating the perturbative evolution factor c' k (jj, pt ) of Eq. (|0|) at NLO in the 
anomalous dimension is well under control. This in turn requires an estimate of the size of the 
NNLO contribution to c' k (fX pt ). To this purpose we have re-computed c^(jU pt ) with two different 

values of the NNLO anomalous dimensions y^ 2 ': in the first case we set \ j k = |){ /){ |> 

1(2) 

in the second case, we guess y k K ' by performing a one-parameter fit to the SSF with sy,S2 fixed by 
perturbation theory, and equating the resulting value of 53 to its perturbative expression 



1(2), n. [ /(0) 1(1) , /(l) , , f(0j 

s 3 = Y k y ; ln2 + Yk Yk + 2b oY k + b iYk 



+ 



(ln2) 2 + 
(ln2)\ 



(6.3) 
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For the operator J2[, we find that in either case the central value of the combination c^(;Uhad) = 
c^/ip^L^jUpt, /ihad) changes by a small fraction of the statistical error, of the order 0.1-0.3 stan- 
dard deviations (depending on the renormalization scheme). For the operator J2' 2 , which carries 
relatively large NLO anomalous dimensions, the effect can be as large as 0.8-1.0 standard devia- 
tions. Therefore, we add to ^(/ihad) a corresponding systematic uncertainty of this order. It has 
to be stressed that the impact of this extra uncertainty at the level of the B-B mixing amplitude is 
not particularly worrying, since the matrix element of £i' 2 enters the latter only at 0(a s ), when the 
static-light theory is matched to QCD. The results for the SSFs and the operator RG running in the 
reference schemes (see the end of Section |3|) are illustrated in Figure ^. 




m/a m/A 

Figure 2: Upper plots: the step scaling functions O12 (discrete points) as obtained non-perturbatively. The 
shaded area is the one sigma band obtained by fitting the points to a polynomial. The dotted (dashed) line 
is the LO (NLO) perturbative result. Lower plots: RG running of =2j 2 obtained non perturbatively (discrete 
points) at specific values of the renormalization scale fi, in units of A. The lines are perturbative results at 
the order shown for the Callan-Symanzik j3 -function and the operator anomalous dimension 7. 



7. Matching to hadronic observables 



The RGI operator is connected to its bare counterpart via the total renormalization factor 
. We stress that Z' kRG1 (go) is a scale-independent quantity, which more- 
over depends upon the renormalization scheme only via cutoff effects. Indeed, it depends on the 



Z (.rgi(So) of Eq. (L6 
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k 


s 


^(jUhad) 




b k 


Ck 


1 


1 


0.777(17) 


0.5731(11) 


-0.171(11) 


0.082(25) 


2 


3 


0.675(12) 


0.7258(14) 


-0.061(14) 


0.016(33) 



Table 1: Running c k (pLhad) (with /Xh a( j = (2L max ) ) and fits to the total renormalization factor of Eq. \.6 



particular lattice regularization chosen, though only via the factor ^'(go^AVd), the computation 
of which is much less expensive than the total RG running factor c' k (lJ.had)- 

We have computed ^'(go^/Xhad), AVd = V(2^max) non-perturbatively at four values of j8 for 
each scheme and four-fermion operator, and for the four different static actions under considera- 
tion. The total renormalization factors are obtained upon multiplying by the corresponding running 
factors on third column of Table |l[ Polynomial interpolations of the form 

ZU>i(so) = a k + h(P - 6) +c k (P - 6) 2 , (7.1) 

can be subsequently used to obtain the total renormalization factor at any value of j8 within the 
covered range (6.0 < /3 < 6.5). We provide in Table [I] the resulting fit coefficients for the HYP2 
action in our reference renormalization schemes. These parametrizations represent our data with 
an accuracy of at least 0.3%. The contribution from the error in the RG running factors of Table |] 
has not been included: since these factors have been computed in the continuum limit, they should 
be added in quadrature after the quantity renormalized with the factor derived from Eq. ( f7Tl| ) has 
been extrapolated itself to the continuum limit. 



8. SF correlators for the bare matrix elements 



In order to simulate the physical matrix elements needed for the B-parameter, we adopt a 
formalism similar to the one described in the previous sections, where the heavy quark field is 
described by the HYP2 static action, while the light quark field is described by the tmQCD action 
including the Sheikoleslami-Wohlert term with non-perturbatively defined c sw . The interpolating 
operators of the B s - and B v -mesons are provided by the boundary sources E s /, and E'r. Correlation 
functions are then constructed by inserting a bilinear or a four-fermion operator in the bulk of the 
SF. Accordingly, the building blocks of the computation are given by 

fx(xo) = — |" L(Xhs(x)L sh ) , f x (xo) = — . ( 8 - D 

Z X Z X 



F Y (x )=a 3 ^(L' s -J hs - hs (x)I. sh ) , (8.2) 



where X = Aq 31 , V stat and Y = J2[ 2 . To be precise, the extraction of the B-parameter of Eq. (1.1) 
requires that the matrix elements of J2[ 2 be normalized by the square of the decay matrix element 
of the By-meson mediated by the static axial current. Since the latter is rotated at full twist into 
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a linear combination of axial and vector currents, we obtain the single contributions to the B- 
parameter from the plateau region of the ratios 



where 



*.-(*<>) = g 



8[2h A _ iV (x )][2h' A _ iV (T-xo)] ' 
1 



1,2 



(8.3) 



r^stat y'Stat/ \ ystaX 

^2 L A.RGL/A \ X 0) -^V.RGU 



BfTM ■ (8-4) 

The RGI axial constant Zff Gl has been non-perturbatively computed in Jl4|, |9|]. The scale indepen- 
dent ratio Z^f 01 /Zff Gl is taken from [15]. We performed simulations at j8 = 6.0,6.1,6.2 with the 
strange quark mass set to physical values as in [16]. Lattice parameters are collected in Table ^[ 

9. Analysis of the excited state contaminations 

The standard way to identify a plateau interval for a three-point correlation function such as 



Eq. (p.2) is to analyse the exponential decay rate of the corresponding meson propagator h A _ iV , 
obtained via the binding energy 



oE eS (xo) 



1, f K_ iW (xp- 

2 \fc A -iv(*0- 



(9.1) 



p 


rxL 3 




K 


n 


6.0 


32 x 16 3 


0.135196 


0.135181 


0.028669 


6.1 


38 x 24 3 


0.135665 


0.135650 


0.028532 


6.2 


44 x 24 3 


0.135795 


0.135785 


0.022890 



Table 2: Lattice parameters 



This procedure may work only provided that the 
lowest value Xq, at which the fundamental state 

is numerically isolated, fulfills the condition x™ m <T/2. Correspondingly, the interval [;e™ n ,r — 
x™ m ] can be certainly used to extract the plateau value of the three-point correlator. Unfortunately, 
this is not the case, as shown in Figs. |3| and ^| (left): due to the small mass gap between the lowest 
and the first excited states in the static-light channel, the plateau starts at about the middle of the 
lattice. Simulations at larger time extensions are increasingly expensive owing to the exponential 
rise of the noise-to-signal ratio related to static propagators. Irrespective of this, the observables /?, 
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Figure 3: The binding energy and the contributions to the Z? s -parameter at j3 = 6.0. 
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are characterized by a very flat time dependence; examples are provided by Figs. || and || (right). In 
order to understand this behaviour, we perform an expansion of Eqs. ( |8.l| , ^ ) through the insertion 
of complete sets of Hamiltonian eigenstates. Assuming that the excited contributions in the vacuum 
channel may be disregarded, one easily arrives at the representation 



R i (x )=B i 



oo g("> m ) 



.(!) 



(n,m)/(0,0) 



fnmSn 



(9.2) 



where A 



(B) 



— is the energy gap between the n-th Hamiltonian eigenstate \n,B) with the 



quantum numbers of a static B s -meson and the fundamental state. Moreover, 



B 



(n,m) 



(n,fi|^-|m,fi) 



fm 



f(«,B|Ao|0,0)(0,0|Ao|m,B) ' 
(i B \n,B}(m,B\i B ) (n,fl|Ao|0,0)(0,0|A |m,B) 



(9.3) 
(9.4) 



(i B \0,B)(0,B\i B ) ' bnm (0,5|A |0,0)(0,0|A |0,B) " 
Here represents the SF boundary state corresponding to the action of the bilinear sources 
Eq. ( [IO| ) on the vacuum. In particular, one should observe that By''™' represents a generaliza- 
tion of the B-parameter, describing the particle-antiparticle mixing of excited states. Numerator 



_. _ R (n,m) /R (0,0) 



and denominator of Eq. (9.2) look quite similar. They only differ by the weighting coefficients 



The hypothetical condition z n m ~ 1 would act on the ratios R{(xo) as an additional damping 
factor of the excited state contaminations, together with the exponential decays due to the mass 
gaps A^. In practice, what really matters for our concern is the first excited contribution, because 
already the second excitation is reasonably expected to compete with the ++ glueball (r^mc ~ 
1.7 GeV). Data suggest that zio could be quite close to one. 
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Figure 4: The binding energy and the contributions to the Z? s -parameter at j3 = 6.2. 

Since the quantum states \B,n) and \B,0) differ only by their mass, we are led to speculate 
about the mass dependence of the generalized B-parameters B " ,m \ Although quantitative state- 
ments are highly non-trivial, it is not difficult to identify at least one extreme situation where the 
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limits z nm = 1 could be realized. This is a scenario in which B) ' is weakly dependent upon the 
mass of the external states. If B^ 1 '" 1 ' is close to fij ' ' , then their ratio will be close to one. This 
picture imposes no restrictions on the value of By ' . 

An apparently different possibility is represented by the VS A, which implies B\ ' = 1, and 
consequently z, m „ = 1. Though speculative, it is not unreasonable that the violation of the VSA 
depends weakly upon the mass of the external states and is responsible in the end for the realization 
of the above-mentioned scenario. 

A quantitative check of the suppression of excited state contaminations in ^1,2(^0) would 
be provided by the level of stability of the observed plateaux under a variation of the boundary 
interpolating operators. In the framework of the Schrodinger functional this possibility is explored 



via the introduction of boundary wave functions like in [ 17]. Unfortunately, the computational price 



required at present for the practical implementation of this smearing technique amounts to giving 



up one of the boundary summations of Eq. (p. 1|), with a corresponding increase of the statistical 



error by a factor of vZA The situation is even worse with a three-point correlator such as Eq. (8.2), 
which has interpolating sources on both boundaries. In this case the introduction of smearing wave 
functions increases the statistical noise by a factor of L 3 , which makes the check useless. This 
problem can be hopefully overcome through the implementation of a SF all-to-all propagator like 
proposed in Jl8|, [l9|, 55]. This is currently under way. 



10. A two-state stochastic model 

In order to have a qualitative view about the impact of large deviations of zio from one on 
the time dependence of Ri(xo), we consider a two-state model. Here, the binding energy and the 
contributions to the B-parameters are described by the stochastic variables 

1 f 1 - oe _A ^ 0_1 ) 1 
£ (e,,oT)=e + -log| 1 -^ TTy | , (10.1) 

f( '* r) ^ i-p[Ui + .^-.)| ' (10 " 2) 

where p and A are differently distributed random coefficients, while z and e parametrize zio and 

(B) 

Eq . Obviously, A is meant to represent the energy gap of the first excited state. From a two-state 
analysis of E e ff(xo), it is roughly known that aA 0.22(3). Therefore, we model this variable 
according to a Gaussian distribution probability, i.e. 

P(A) = — ^=exp , (A,(T A ) = (0.22,0.03) . (10.3) 

On the other hand, p is supposed to represent the product of the matrix elements /10 and gio, 



defined in Eq. (3.4). Choosing a distribution probability for this variable is delicate, because we are 
largely ignorant about the projection of the SF boundary state |i#) onto the first excited state 1 1,2?) 
and the decay constant of the latter. We can heuristically expect that 

(0,S|A |0,0) m B f B 
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Nevertheless, if we believe that is well projected onto |0,B), then /io ~ 0. In this case we 
should choose a probability distribution of p peaked around p = 0. By contrast, if we believe that 
is a balanced mixture of |0,B) and it follows that /io ~ 1. It makes sense to assume 

a given sign for p and not to allow for fluctuations of the opposite sign. A flexible distribution 
probability allowing for a definite sign is the Log-normal distribution, defined by 

P(p;p,o p ) = 7^^P{- ' 2 \ • (10-5) 

Having produced N samples {A,}, = i ^ and {pi}j=i...N of A and p, we approximate the ensemble 
averages of h(z,xq,T) and e(e,xo,T) via 

S{e^T) = (e(e, X0 ,T)) ^ e + — £log| j , (10.6) 

1 W 1 _ r e -A,(.t ) _|_ e -A i (7'-.to)] 

*(z,*,r) = {P(Z , X0 ,T)) ^I ^^^^.^ • (10-7) 

The ensemble averages <f and ^ are now functions of the distribution parameters p and a p , which 
can be varied in order to change the shape of the distribution. One of the worst cases we considered 
is the one corresponding to (p = —l/8,o p = 1/4), i.e a Log-normal distribution peaked around 
exp(— p) ~ 0.88. In the spirit of the two-state model, such a distribution describes a large overlap 
between the interpolating boundary state and the first excited one. As shown in Fig. BL the binding 
energy resembles very closely the one of Fig. ||. The shapes of the B-parameter corresponding of 
different choices of z suggest that zio could be very close to one in the real case, thus supporting 
our interpretation in terms of the VSA. 

11. Conclusions 

B° — B° mixing remains among the most important processes that are required to pin down 
the elements of the CKM matrix precisely. However, in order to constrain the unitarity triangle 
sufficiently well and to look for signs of new physics, theoretical uncertainties associated with 
hadronic effects must be reduced. In this talk we have reported on a new strategy for the com- 
putation of the heavy-light B-parameters in lattice QCD, based on tmQCD and HQET. Its main 
advantage is the exact absence of mixing under renormalization, which plagues standard Wilson 
fermions, at the same computational cost of a Wilson-type regularization. We have described our 
fully non-perturbative calculation of the relations between parity-odd, static-light four-quark op- 
erators in quenched lattice QCD and their renormalized counterparts. We have also described our 
first experiences with the computation of the bare matrix elements for the B s — B s mixing and their 
excited state contaminations. Before attempting a continuum extrapolation of the matrix elements, 
a deeper analysis of the excited state contributions has to be done and we hope that the all-to-all 



propagator, like proposed in [ ]18| , \l% [20J will be of great help there. 
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